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In-medium effects of transverse-mass distributions of quarks and gluons are con- 

bb" 

sidered assuming a possible local equilibrium for colorless quark objects like mesons 

o" 

CN ■ and baryons created in central A+A collisions. It is shown that the average square 

qj , of the transverse momentum for these partons grows and then saturates when the 



(N 
(N 



initial energy increases. Within the quark-gluon string model it leads to the energy 
dependence of hadron transverse mass spectra which is similar to that observed in 
heavy ion collisions. Comparison with other scenarios is given. 

§ ; 1. INTRODUCTION 

Searching for a new physics in heavy-ion collisions at AGS, SPS and RHIC energies has 

led to intense theoretical and experimental activities in this field of research [1]. In this 

CO ■ respect the search for signals of a possible transition of hadrons into the QCD predicted 

phase of deconfined quarks and gluons, quark-gluon plasma (QGP), is of particular interest. 

Among many signals of the QGP formation, one of the earliest is based on the relation of 

the thermodynamic variables, temperature and entropy, to the observable quantities, the 
• i-^ . 

average transverse momentum and multiplicity, respectively. Shuryak [2, 3] and van Hove [4] 

hi , 

were first to put forth arguments that the transverse momentum as a function of multiplicity 
reflects general properties of the equation of state of hot and dense matter. In the case of the 
first order phase transition this relation is quite specific due to inevitable formation of the 
coexistence phase of quark-gluon plasma and hadrons, the mixed phase (MP). Experimental 
detection of the QGP phase and the MP in A+A collisions is a non-trivial task because of 
smallness of the space-time volume of the hot and dense system and possible contributions 
of hadronic processes simulating signals of QGP and MP formation. Nevertheless, recent 
experimental study of the transverse-mass spectra of kaons from central Au+Au and Pb+Pb 
collisions revealed "anomalous" dependence on the incident energy. The effective transverse 
temperature (the inverse slope-parameter of the transverse mass distribution at the mid- 
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rapidity) rather fast increases with incident energy in the AGS domain [5], then saturates at 
the SPS energies [6] and approaching the RHIC energy region [7] increases again. In agree- 
ment with expectations this saturation was assumed to be associated with the deconfinement 
phase transition and indication on the MP [8, 9]. This assumption on the first order phase 
transition was indirectly confirmed by the fact that two independent transport models, based 
on hadron-string dynamics, UrQMD and HSD, failed to reproduce the observed behavior of 
the effective kaon temperature [10, 11]. Though the calculated transverse mass spectra are in 
reasonable agreement with the experimental results for pp, pA and light nuclei (central C+C 
and Si+Si) collisions, both models strongly underestimate the transverse energy T above en- 
ergy 5 AGeV and are not able to describe the rapid increase of the inverse slope parameter 
for collisions in the AGS domain and its subsequent flattening in the low energy SPS regime. 
The failure was attributed to a lack of pressure. This additional pressure must be generated 
in the early nonhadronic phase of the collision because the strong hadronic interactions in 
the later stages do not produce it [10, 11]. The anomalous effective-temperature behavior 
was quite successfully reproduced within hydrodynamic model with equation of state involv- 
ing the phase transition [12]. However, this result is not very convincing since to fit data 
the required incident-energy dependence of the freeze-out temperature should closely repeat 
the shape of the corresponding effective kaon temperature and thereby the problem of the 
observed anomalous inverse-slope dependence is readdressed to the problem of the freeze-out 
temperature. 

In addition, as was shown in [13], these experimental data within the considered energy 
range up to ~ 160 AGeV may be reasonably described in the 3-fluid relativistic hydrody- 
namic model with the pure hadronic equation of state. Note that this 3-fluid model also well 
describes a large variety of global observables [14]. This implies that a hydrodynamic collec- 
tive motion in the expansion stage, which is absent in transport codes, plays an important 
role in the anomaly discussed. It is of interest that recently a step to taking into account 
some collective motion in the microscopic description has been done in [15] where a hadronic 
transport model has been extended to include three-body collisions. The agreement between 
the transport model and experiment for the kaon transverse-temperature excitation function 
was noticeably improved. 

In this paper we propose another way to introduce a collectivity effect in nuclear system 
via some in-medium effect. We consider the temperature dependence of quark distribution 
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functions inside a colorless quark-antiquark or quark-diquark system (like meson or baryon, 
h) created in central A+A collisions. A contribution of this effect into transverse momentum 
spectra of hadron is estimated and it is shown that it results in larger values of the inverse 
slope parameter and therefore in broadening of the transverse mass spectra. 



2. QUARK IN A HADRON EMBEDDED IN EQUILIBRATED MATTER 

Let us assume the local equilibrium in a fireball of hadrons whose distribution function 
can be presented in the following relativistic invariant form: 

f£ = C T {\ ± exp((p h -u-fi h )/T)}- 1 , (1) 

where ph is the four-momentum of the hadron, the four- velocity of the fireball in the proper 
system is u = (1, 0, 0, 0), the sign "+" is for fermions and " — " is for bosons, \ih is the baryon 
chemical potential of the hadron h, T is the local temperature, and Ct is the T-dependent 
normalization factor. The distribution function of constituent quarks inside h which is 
in local thermodynamic equilibrium with the surrounding nuclear matter, f^(x, p t ), can 
be calculated using the procedure suggested for a free hadron in Ref.[16]. So, the quark 
distribution in a free hadron moving with the momentum P is given by the convolution 

ft(Pz, Pt) = J dp lz d 2 p lt q v (p z , pt) q r (pu, Pit) S(p z + p lz - P) 5 (2) (pit + Pt) , (2) 

where q v (Pz,Pt) is related with the probability to find the valence quark with longitudinal 
momentum p z and transverse momentum p t in the hadron, whereas q r (pi z ,Pit) is the prob- 
ability that all the other hadron constituents (one or two valence quarks plus any number 
of quark-antiquark qq pairs and gluons) carry the total longitudinal momentum p\ z and the 
total transverse momentum pu- 

Introducing the new variables x = p z /P,Xi = p\ z /P, which at large P become the 
relativistic invariant light cone variables, we can rewrite eq.(2) (for example, for the valence 
quark distribution) in the following equivalent form [17, 18]: 



fg v (x,Pt) = J Q dxi J d2pit ^( X 'P*) 9r(a:i,Pit) 5(x + X! - 1) 5 (2 \p lt + p t ) , 



(3) 



Assuming the factorization hypothesis q v ,r(x,Pt) = <lv,r(x) g v ,r(pt) and integrating eq.(3) 
with respect to dx\ and d 2 p u we have 

/£(*,ft) = /£0O$t(pt), (4) 
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where 

ft(x) = Qv(x) q r (l - x) , gg v (p t ) = 9v(Pt) 9r(-Pt) (5) 
For simplicity one can choose the normalized ^-distribution g^(p t ) in the Gaussian form 



7<? 

Then we get 



g h M) = - ex P(-w 2 ) • (6) 

71 



!h(jH) = U-Vt) = (^) 1/2 exp(-|p t 2 ). (7) 

Similarly to eq.(2), one can present the quark distribution inside the hadron h to be in 
local thermodynamic equilibrium with the surrounding matter, f^(x, pt), in the convolution 
form as 

fq(p*,Pt) = J dp lz dp hz d 2 p lt d 2 p ht q v (Pz,Pt) 9r(pi*,Pit) fh(Pht,Pht) (8) 

x 5(p z + p lz - p hz ) 5 {2) (p t + pu -Pht) ■ 

In contrast to eq.(2) the distribution of the hadron h in a fireball f£ is included in 
eq.(8), therefore we integrate over the longitudinal and transverse momenta of h. In 
the general case, the functions q v (piz,Pit) and q r (piz,Pit) entering into eq.(2) may dif- 
fer from q v (piz,Pu) and q r {Piz,Pit) in eq.(8). Let us introduce the Feynman variables 
x = 2p*/\fs i , xi = 2p\ z /\fs i ) %h = 2p* hz /\fs i , where P* z ,p\ z ,Ph z are the longitudinal mo- 
menta and s' is some characteristic energy squared scale. Then eq.(8) can be reduced to the 
following form similar to eq.(3): 

ftv( x >Pt) = / dx i / dx h / d 2 p lt d 2 p h t q v (x,p t ) q r (xi,Pit) (9) 
Jo Jo J 

xfh(xh,Pht) dix + xx-Xh) 5 {2) (p lt + p t - Pht) ■ 

If the factorization hypothesis q v>r (x, Pt) = q v ,r(x) g v ,r{pt) and the Gaussian form for g v , r (pt) 
are assumed, 

g x Apt) = (^y 2 ex P (-|p t 2 ), (10) 
we can get the following expression for f^(x,p t ) (see APPENDIX) : 

ff(x,p t ;T) = 2m h TC T exp( h fc )(lH )- / dx h q v (x) q r (x h - x) 

1 171^ TT J x 

x f q (x h ) exp(-f q (x h )p 2 ) , (11) 
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where rrih(xh) = \/ m h + Xhs'/4, by assumption 7 g = 7 g and 

f / n _ 7g(l + 7^K)^/2) n v 

l + lq m h {x h )T " 

From the normalization relation for f^(x, pt] T) 

dx j d 2 Pt ff(x,p t ;T) = 1 (13) 
we can find that C T = (2m h T(l + T/m h ))~ l exp[(m h — fj, h )/T], q v (x) = q v (x) and 

pi pl—X 

q r (l -x) = dx h q r (x h - x) = q r {y)dy . (14) 

Jx JO 

Then eq.(ll) reads 

1 r l ~ x _ 
fq{x,p t ]T) = - dxi q v (x) q r {xi) T q (xi + x) exp(-T q (xi + x)p 2 ) . (15) 

^ Jo 

As is evident from eq.(15), at vanishing temperature T = the quark distribution 
reproduces the quark distribution in a free hadron. 

Using the quark distribution f^(x,p t ) we can calculate the average value for the transverse 
momentum squared for the quark < p 2 t >^ in h as a function of x and s^n- 



, , v A _ J f q A (x,Pt)p 2 t d 2 Pt 
<P ^ X)>h ~ f f q A (x,p t )d*p t 



(16) 



Applying eq.(15), we have 



<P,»>„=C j / o ( 17) 

where < p\ > q = l/7 g is the average transverse momentum squared for a quark in the free 
hadron and 

-l-x 



C' 1 = / dxxq^xx) . (18) 
Jo 

At x ~ eq.(17) can be presented in the following equivalent form: 



2 , . ..... i /, /'' , ^■(•'•|)(< P? >J • ( v '"'/; • IV) 



<pi t (x~0) = C / rfzi /v 9 vv — -, (19) 

9,41 ^ * A 1 l + ( v /m 2 h + xfs'/A)T/(2<p 2 >^ { 1 

where 

C^ 1 = / dxxqrixx) . (20) 
Jo 
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In the general case the quark distribution in a nucleon fq{x) at low momentum transfer 
when its Q 2 QCD evolution can be neglected is presented as follows: 

f q N (x) = C q x a (l-x) b , (21) 

where C q is the normalization factor. The parameters a and b can be extracted from the 
deep inelastic scattering or calculated within some quark models. Therefore, according to 
eq.(14), the function q r (xi) entering into eqs. (17,18) has the following form: 

q r ( Xl ) = bxl' 1 . (22) 

Usually, 6 > 1.5 (see below). So q r (xi) falls down very fast when x\ decreases from 1 to 
and may be taken out of the integral in eq.(19) at x 1 — 1 : 

2 A <Pt >q + T V m h + S 7 4 

< P s>t (* * 0) > Kanr , ~ - - T ^ mi = = = 4/(2 - - — , (23) 



As is evident from eq.(23) the quantity < p qt >^ depends on the energy V s' and temperature 
T. At any nonzero values of T it grows when ^fs* increases (because its derivative with respect 
to is positive) and then saturates, reaching the asymptotic value about 2 < p 2 >^ at 
high energies. The integration over dx\ in eq.(19) does not change this result qualitatively. 

To calculate < Pg t (x) more accurately, we have to know the functions q v (x) = q v (x) 
and q r (l — x) which is related to q r (y), see eqs. (14), (22). Parameters of these distributions 
can be defined by application of the quark-gluon string model (QGSM) [21, 22] based on the 
Regge asymptotic of quark distributions in a nucleon and 1/N expansion in QCD [19, 20] 
(N is the number of flavors or colors). According to this model, the nucleon consists of a 
quark, diquark, and quark-antiquark see (qq). For example, for the valence w-quark in the 
proton [22] we have 

a = -a R (0) ; b = a R (0) - 2a N {0) (24) 

where q;_r(0) = 1/2 is the reggeon intercept of the Regge trajectory, Q;Ar(0) = —0.5 is the 
intercept of the nucleon Regge trajectory. Comparing eq.(21) and eq.(5) one can find that 

Uv (x) = x- aR{0) ; u r (l -x) = (l- x ) a n(°)"^N(o) ( 25 ) 

The similar form can be obtained for the ^-distribution of the valence o?-quark in the nucleon. 
Then, using eq.(19), the average transverse momentum squared can be estimated for w-quark 
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inside the proton which is in local equilibrium in the fireball < p 2 ut This quantity at 
x ~ is presented in Fig.l as a function of \fs i . The behavior of this quantity given by 
eq.(19) is similar to that discussed above in respect to the approximate expression (23). 

3. TRANSVERSE MOMENTUM SPECTRA OF HADRONS FROM CENTRAL A + A 

COLLISIONS 

Let us now estimate the p t distribution of the hadron hi produced after collision of two 
hadrons one of them is locally equilibrated in a fireball. We shall explore the QGSM based 
on 1/N expansion in QCD. Actually, this is the expansion of the QCD amplitude in different 
topologies [19, 20]. 1 The first order term is the so-called planner graphs corresponding to 
the one-Reggeon exchange diagrams in the t-channel of the hadronic process. The second 
order term of this expansion is the so-called cylinder graphs related to the one-Pomeron 
exchange diagrams. The last ones make the main contribution to inclusive spectra of particles 
produced in inelastic hadronic processes. In Fig. 2 the cylinder graphs for inelastic meson- 
nucleon (left diagram) and nucleon-nucleon (right diagram) inelastic processes are presented, 
see also Ref.[20]. According to this model, the colorless strings are formed between the 
antiquark/quark (q/q) in the colliding meson and the quark/diquak (q/qq) in the colliding 
nucleon (left diagram of Fig. 2), then, after their break, qq pairs are created and fragmentate 
into a hadron h±. The contribution of the cylinder graph (right diagram of Fig. 2) to the 
inclusive spectrum = J^ h of the hadron h\ from the collision of two nucleons can 
be presented as follows: [22, 26]: 

pt(x,p t ;T) = a^^x^pt-T^ix-^t-^/F^pt-T) + (26) 
F^(x + ,p t -T)F^(x-,p t -T)/F q (0,p t -T)} , 

where <J\ is the cross section of the 2-chain production, corresponding to the s-channel 
discontinuity of the cylinder (one-Pomeron) graph. It is usually calculated within the quasi- 

1 This model, the Quark-gluon String Model (QGSM) [21, 22] or the Dual Parton Model (DPM) [23] differs 
from the Lund String Model, see for example Ref.[24]. 
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eikonal approximation [25]; 

F^(x ± ,p t ,T) = f ^ I d2 Vu d 2 p 2t $ w) (*i,p u ;T) G^ ?) (^,p 2t ) 

r, J X+ J X l 

flavors ± 

x 5^{p lt + p 2t -p t ) • (27) 

Here x± = \ ( \/ x\ + x 2 ± x) and x t = 2m hlt / ^s hhl where s' has been associated with the 
energy squared s hh of colliding hadrons, m hl t = sjm 2 hi + p\ t ; z^G 11 ^ = D q ^ qq) is the frag- 
mentation function (FF) of the quark q (diquark qq) into the hadron h\ . Actually, the inter- 
action function F q ^(x + ,p t , T) corresponds to the fragmentation of the upper quark/diquark 
( see right diagram in Fig. 2) into the hadron h±, whereas F^^x^jPt, T) corresponds to the 
fragmentation of the down diquark/quark into h\. We can get a similar expression for the 
Pt-spectrum of the hadron h\ produced in the meson-nucleon collision by replacing the di- 
quark qq with the antiquark q. To calculate the transverse momentum spectrum of the 
hadron h\ in the central rapidity region, eq.(27), one needs to know the ^-dependence of the 
fragmentation function D^ 1 . We assume the same Gaussian dependence as eq.(6). However, 
the slope of this p t dependence 7c can differ from the slope 7(J for the constituent quark p t 
distribution. For the mid rapidity region, x ~ 0, eq.(26) can be rewritten in the following 
form: 

pl(x ~ 0,p t ;T) = a^F^x ~ 0,p t ;T)) + F%((x ~ 0,p t ;T)} , (28) 

We have the following approximate equation for the averaged transverse momentum 
squared of the hadron hi produced in the interaction of the nucleon with the nucleon locally 
equilibrated in a fireball: 

2 aa F 9 (a:~l) + 7 C 11 

< Ph x t > NN,appr— = fTl 77 + 77" ' (^9) 

T q (x ~ l) 7c T q (x ~ 1) 7c 



By substitution oit q (x ~ 1) given by eq.(12) at x^ — 1 we finally have 

< p\ > N q +T^m% + s hh /A 
l + T v /m5 f + % /4/(2<p i 2 >f) 



o * A „ < Pt >l +T y/rnj + %/4 < pj > N q 

^ Phxt > NN,appr — 1 , „ / — s — i 77 / /o 2 N\ r- ' ^ ' 
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where r = 7c / 7g . 

From comparison of eq.(30) and eq.(23) one can see that the energy dependence of the 
transverse momentum squared for the produced hadron, < p\ lt >N J N,appr.i i s qualitatively 
similar to that of the quark inside the hadron which is in local equilibrium in the fireball, 
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< p 2 qt > A . At large values of 7 C >> j q , the second term in eq.(30) can be neglected and 

< Ph u t >J NN, apP r. as a function of y/shK and T is close to < p 2 t > A ^ appr _ given by eq.(23). 
In this case the hadron spectrum copies the quark spectrum. In fact, 7 C is larger than ^ q 
by factor r = 3 — 4 as follows from inclusive jvspectra of hadrons produced in hadronic 
processes analyzed within the QGSM Ref.[26]. 



4. RESULTS AND DISCUSSION 

We estimated the average value of transverse momentum squared for K + -mesons produced 
in nucleon-nucleon < p 2 K+ t > A t A j an ^ pion+nucleon < p 2 K+ t > A £ interactions of two hadrons 
one of them is thermodynamically locally equilibrated in a fireball created in the central 
A + A collision. This quantity was estimated as a function of y/Shh at T = 150 MeV. for two 
cases when r y c » lq and 7 C = ?>^ q [26] using eq.(29) for the inclusive spectrum at x ~ 0. It 
is presented in Fig. 3. 

As is seen from Fig. 3, the results obtained are sensitive to the mass value of the hadron 
which is in local equilibrium with the surrounding nuclear matter at y/Shh < 10 GeV. 

Collectivity due to thermal effects results in growth of < p 2 K+ t > AA with energy for 
colliding hadrons in central A-A collisions similar to that for unbound quarks in hot mat- 
ter, < p 2 t > A . At larger energies yjshh this quantity saturates. The saturation value for 
this quantity depends on the hadronization mechanism of quarks/diquarks to hadrons. A 
simple exponential estimate of ^-spectra for produced hadrons hi is used to parameterize 
experimental data: 



dN m hlt 



2— 7-|j/=o = C exp{--^) , (31) 



where the parameter T* is extracted from fitting experimental data. There are data on the T* 
values for different hadrons and different m t domains: "lowp t " when m hlt — m hl < 0.6 GeV, 
and "high p t " , 0.6 < m hlt — m hl < 1.6 GeV (see for example Refs.[5-7]). At low p t the m t 
spectrum given by eq.(31) can be presented in the following approximate form: 



^^-1^=0 ^ Cexp(-m hl )exp(- Vhx \ _ ) (32) 



dmfdy ' IrrihiT* 

Actually, at small transverse momenta 2m^T* ~< p\ ut > AA , where < p\ l t > AA is the 
transverse momentum squared for the hadron h\ produced in central A + A collisions. For 
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fT-mesons < p 2 Kt > Aj4 ~ T** GeV/c 2 . The experimental data on the inverse slope T* for 
f^ + -mesons produced in central Au+Au (Pb+Pb) collisions as a function of the incident 
energy per nucleon show [5, 6] that T* grows and saturates later on. 

Our results presented in Fig. 3 qualitatively demonstrate similar behavior for < p 2 K+ t 
as a function of y/Shh- One should emphasize that here y/Shh is the energy of a pair of 
colliding hadrons and it is not related directly to the initial energy of colliding heavy ions. 
These results are only an illustration of the collective effects assumed. In a real case, such 
binary interactions occur between various hadrons in a large range of temperatures. 

As was noted in the introductional part, the broadening effect for mi-spectra of the 
hadron produced in central A+A collisions observed at AGS, SPS and RHIC energies can 
qualitatively be explained using the assumption of possible creation of the QGP and a 
co-existing phase of quarks and hadrons [8, 9]. From available hydrodynamic calculations 
of the transverse inverse-slope excitation function [12, 13] it is not yet clear whether this 
behavior can be considered as a signal of the phase transition into the QGP. Microscopic 
transport models [10, 11] taking into account formation and decay of strings as well as the 
multiple rescatterings of hadrons are definitely not able to describe these data. An attempt 
to enhance the rescatterings was undertaken in Ref. [28] by including the Cronin effect. In 
the transport approach enhancement of the intrinsic quark transverse momentum spread 
< p\ >g is simulated by increasing the average transverse momentum of quarks < p 2 > q 
with the number of previous collisions of primary nucleons N prev as 

< Pt >f=< Pt > q (1 + ctN prev ) , (33) 

where the parameter a ~ 04. Now the description of spectra becomes rather good at the 
RHIC energies, improves essentially at the SPS energy of 160 AGeV, but does not show 
any significant change at 11 AGeV [28]. Consequently, the "pre-hadronic" Cronin effect, 
realized via eq.(33), is not responsible for the anomalous behavior of kaon slopes around 
AGS energies. 

Another scenario of collectivity, color rope formation, was proposed in [29]. This color 
rope model assumes that in central A+A collisions several strings are produced, some of 
them on top of each other. The common chromo-electric field created by overlapping K 
single guark/antiquark sources may form a if- fold rope. The spread of hadron transverse 
mass distributions resulting from color string-rope decay is defined by the surface tension 
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parameter k k which is VW times larger than the appropriate parameter for the decay of a 
single string in N — N collisions, k k = VWk. In a certain sense, this scenario is opposite to 
the proposed one: changing originates not from the initial "pre-hadronic" level, but rather 
comes from the final state as in-medium modification of the string break function. Similar 
effect for overlapping strings was estimated in Ref.[28]. Only a small increase in the inverse 
slope parameter at AGS energies was found because the string densities are low. At SPS 
and RHIC energies the model gives hardening of the spectra by about 15% [28]. 

5. CONCLUSION 

We have found that the quark distribution in a hadron depends on the fireball temperature 
T. At any T the average transverse momentum squared of a quark grows and then saturates 
when \fs i increases. Numerically this saturation property depends on T . The modification of 
of initial quark distributions leads to a similar energy dependence for the average transverse 
momentum squared < p\ l t of the hadron h\. The saturation property for < p\ i t 
depends also on the temperature T and it is very sensitive to the dynamics of hadronization. 
As an example, we estimated the energy dependence of the inverse slope of transverse mass 
spectrum of i^-mesons produced in the interaction of two hadrons in the fireball created in 
central A+A collisions. It is qualitatively similar to the incident energy dependence of this 
quantity observed experimentally. We guess that our assumption on the thermodynamical 
equilibrium of hadrons given by eq.(l) can be applied for heavy nuclei only and not for the 
early interaction stage. 

From the above discussion we see that the observed anomalous behavior of the kaon inverse 
slope in central A+A collisions is still puzzling. There are several scenarios which can be 
valid in various degrees, however, the final consistent solution of this puzzle is still absent. Its 
solution can partially be due to the proposed thermal mechanism of collectivity of hadrons. 
For final decision this and other effects should be taken into in a dynamical transport model. 
The solution of this "step-like" puzzle is an important point in the scientific programs on 
the future heavy-ion accelerators FAIR [30] and NIC A [31] 
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Figure Captions 

Fig.l. The energy dependence of the average transverse momentum squared for the w-quark 
in a proton in nuclear matter at temperature T. 

Fig. 2. The cylinder graph for the inelastic meson-nucleon processes (left) and the cylinder 
graph for the inelastic nucleon-nucleon reaction (right) [20] 

Fig. 3. The average transverse momentum squared of the lC~-meson produced from the 
interaction of two hadrons one of them is in the equilibrated fireball as a function of its 
energy y/Shh at T = 0.15 GeV. Curves 1 and 2 correspond to < p 2 K+ t >jj% and < p 2 K+ t >£fi 
respectively when 7 C >> 7 9 , whereas curves 3 and 4 correspond to the same quantities when 
7c = 37g. Line 5 corresponds to the average transverse momentum squared of K + produced 
in free p + p collisions < p 2 = 0.14 GeV/c 2 . 




Figure 3. 
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APPENDIX 

Let us get eq.(ll) using eqs.(9,7) and eq.(6). Eq.(9) can be also written in the following 
equivalent form 

fg(x, p t , T) = ^-C T exp(j^) J dx h q v (x)q r (x h - x) exp(-T^) x (34) 



X (2^ / exp (~ 7 ^^ X ex P(~ ^ } 6XP( ^' ( * u + ^ 2t " Pt))d 2 k lt d 2 k 2t d 2 b , 



where 7^ = %/2, rrVfJjCh) = m h + x\s' /A. The first integral in eq.(34) reads 

■h = /e.pC- ^y^W g- - 2- reM- Vkl+ f h>X ^) x (35) 



xJo(bk 2t )k2tdk 2 t = 27r ^ 2 + 52)3/2 ^ + wifc(zh) \ja 2 T + b 2 ) exp(-m ?i y / 'a|, + 6 2 ) , 

where ay = 1/T, Jo(bk2t) is the Bessel function of order depending on 6/c 2 t . 

Eq.(35) in the central region when 6 2 < 1/T 2 can be presented in the following approximate 

form: 

Ji = 27rT(m ft (x fc ) + T) exp{-b 2 rh h {x h )T/2) . (36) 

Using now the form for J\ given by eq.(36) we can calculate the second integral in eq.(34) 
f b 2 b 2 m h T ^ .-. 2 ±irj q j qV 2 

Including all the terms staying in front of eq.(35) we get the form for f^(x, p*;T) given by 
eq.(ll). Note that by getting eq.(ll) the term 

{m h (x h ) + T) exp( — ) 

was moved out the integral in eq.(34) at Xh ~ because the exponential function falls down 
very fast when Xh increases from up to 1. 



